Intrinsic spin fluctuations reveal the dynamical response function of holes coupled to 

nuclear spin baths in (In,Ga)As quantum dots 
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The problem of how single "central" spins interact with a nuclear spin bath is essential for un- 
derstanding decoherence and relaxation in many quantum systems, yet is highly nontrivial owing to 
the many-body couplings involved. Different models yield widely varying timescales and dynamical 
responses (exponential, power-law, Gaussian, etc). Here we detect the small random fluctuations of 
central spins in thermal equilibrium (holes in singly-charged (In,Ga)As quantum dots) to reveal the 
timescales and functional form of bath-induced spin relaxation. This spin noise indicates long (400 
ns) spin correlation times at zero magnetic field, that increase to ~5 /.is as hole-nuclear coupling 
is suppressed with small (100 G) applied fields. Concomitantly, the noise lineshape evolves from 
Lorentzian to power-law, indicating a crossover from exponential to inverse-log dynamics. 
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Single electron or hole "central" spins confined in III-V 
semiconductor quantum dots (QDs) are promising can- 
didates for solid-state qubits Although confinement 
suppresses momentum-dependent spin relaxation path- 
ways, it enhances the hyperfine coupling between the cen- 
tral spin and the dense spin bath of ~10^ lattice nuclei 
comprising the QD. These hyperfine interactions dom- 
inate decoherence and spin relaxation at low tempera- 
tures. Within a QD ensemble, each central spin feels 
a different effective (Overhauser) magnetic field from 
nuclei, Bn. Trivially, this leads to rapid nanosecond- 
timescale ensemble dephasing of an initially-oriented en- 
semble of central spins 0]. On longer timescales, how- 
ever, decoherence and relaxation of the central spin 
within each QD occurs because Bn evolves slowly in time 
[ll-[lot. In large applied magnetic fields B, the huge dif- 
ference between electronic and nuclear Zeeman energies 
suppresses flip-flop interactions between the two species, 
and B„ evolves primarily via weak dipolar coupling be- 
tween nuclei. As B— J-O, however, central spins can facil- 
itate mutual "co-flips" between distant nuclei (a process 
involving virtual flips of the central spin) , which changes 
Bn more rapidly and accelerates spin relaxation [l|-[lOj. 

It is precisely in this low-field, intimately-coupled 
regime where the decoherence of the central spins be- 
comes exceedingly difficult to model theoretically. Dis- 
tinctly different timescales and a wide range of dynam- 
ical response functions (exponential, Gaussian, power- 
law) have been postulated, with exact solutions derived 
only under certain limiting assumptions, such as po- 
larized nuclei. Low-field numerical models with un- 
polarizcd nuclei suggest interesting non-exponential dy- 
namics with slow l/log(i) decays highlighting the 
non-Markovian and strongly-correlated evolution of this 
many-body quantum system. 

While groundbreaking experimental QD studies fo- 
cused on electron central spins, considerable attention 



has recently shifted to holes [llHl9j , whose p-type wave- 
functions avoid strong Fermi-contact hyperfine coupling 
to the lattice nuclei. Instead, hole-nuclear coupling oc- 
curs primarily via weaker dipolar (anisotropic hyperfine) 
interactions, reducing Bn by one order of magnitude 
[2l| . Optical studies of QD holes based on re- 
peated initialization [HI . or continuous pumping 
113, have revealed long spin relaxation and coherence 
times in large B. However, studies in the B— >-0 limit, 
where hyperfine interactions are manifest most strongly, 
have received comparatively little attention 17[. More- 



over, the underlying functional form of the dynamical 
response of holes in a spin bath has not been explored. 

As an alternative to conventional pump-probe tech- 
niques, the fluctuation-dissipation theorem suggests an- 
other route to reveal the dynamical response function 
of holes, that is based on passively detecting the spec- 
trum of intrinsic random spin fluctuations of holes in 
thermal equilibrium (i.e., without optical pumping or ini- 
tialization) . This "spin noise spectroscopy" has origins in 
atomic physics 22|, |23| and nuclear magnetic resonance 



2J], and was subsequeritly demonstrated for electrons 
in bulk semiconductors 23, 26|. As recently applied to 
QDs [13], spin noise revealed the precession and rapid 
ensemble dephasing of holes in large B, due to an inho- 
mogeneous distribution of hole g-factors. However, these 
noise-based methods have never been used to detect the 
intrinsic dynamics of central spins interacting with a nu- 
clear spin bath. 

We therefore use a passive optical technique based on 
Faraday rotation to detect intrinsic hole spin fluctuations 
in (In,Ga)As QDs as B—>-0, where coupling to the nu- 
clear spin bath is most important. Crucially, because 
spin correlations are revealed in the spectral domain, this 
approach is well-suited to determine slow dynamical re- 
sponse functions with an accuracy sufficient to achieve 
a novel understanding of coupled spin-bath systems. In 
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FIG. 1. (a) Experimental schematic: The random spin fluc- 
tuations SS'^{t) of resident holes in (In,Ga)As/GaAs QDs 
impart Faraday rotation fluctuations 56 pit) on a linearly- 
polarized probe laser. The power spectral density of this 
"spin noise" is measured with a balanced photodiode bridge 
and digital spectrum analyzer (see also [231). (b) Typical 
spin noise power spectrum of the resident holes at low tem- 
perature (5 K) and zero applied magnetic field [B=Q). The 
400 kHz half- width of the spin noise indicates a long ~400 ns 
correlation time rj* of the hole spins, (c) The same spec- 
trum on a log-log scale. The noise lineshape closely follows a 
Lorentzian, indicating exponentially-decaying hole_S£in cor- 
relations at 5=0, in contrast with recent theories 



contrast with theoretical predictions exponential 
dynamics with long (400 ns) correlation timescales are 
found at B=Q. Using small (100 G) applied fields to sup- 
press a dominant hole-nuclear interaction channel, even 
longer timescales of order 5 /iS are revealed. Concomi- 
tantly, the fluctuation spectrum evolves from Lorentzian 
to power-law, indicating a crossover from exponential to 
inverse-log spin relaxation. 

Figure 1(a) summarizes the spin noise experiment. A 
linearly-polarized probe laser is focused through an en- 
semble of weakly p-type (In,Ga)As/GaAs QDs, where 
~10% of the QDs contain a single hole. Stochastic fluc- 
tuations of the ensemble hole spin projection along the 
sample normal z, SS!^{t), impart Faraday rotation fluc- 
tuations 66p(t) on the probe laser via the usual optical 
selection rules for positively-charged trions. The power 
spectral density of 59p{t) is measured, revealing hole spin 
fluctuations in the frequency domain. Via the Wiener- 
Khinchin theorem, this is equivalent to the Fourier trans- 
form of the hole spin correlation function (S'^(i)5''(0)). 
We emphasize that this passive technique probes the in- 
trinsic fluctuations of the resident hole spins while in 
thermal equilibrium with the nuclear spin bath - to lead- 
ing order the detected spins are not optically pumped or 



excited [28|, in marked contrast to conventional pump- 
probe methods. Moreover, since {S^{t)) is always zero, 
only the intrinsic long-timescale dynamics of hole spins 
are detected in the zero- field limit. In contrast with prior 
work 27|, two crucial differences are implemented here: 
i) low-noise stabilized probe lasers that permit accurate 
and quantitative recovery of the small low-frequency spin 
noise signals that exist in the zero-field limit, and ii) very 
low probe powers (< 0.5 mW) and very large spot sizes 
(^100 fim) to ensure that we operate in a regime where 
the probe laser itself does not influence the measured spin 
noise signals and lineshapes (for details, see [28[). 

Using this approach, we find that the spectral density 
of this hole spin noise in zero applied field consists of a 
single, well-defined peak centered at zero frequency (Fig. 
1(b)). Considerable information is encoded within this 
noise peak: Its half-width F reveals the characteristic 
timescale of hole spin correlations {S^{t)S'^{0)) , and 
most importantly its detailed lineshape directly reveals 
the functional form of the central (hole) spin decay - 
a parameter of considerable theoretical interest |ll4l0|. 
On a log- log scale (Fig. 1(c)), we find this hole spin 
noise closely follows a Lorentzian lineshape over three 
orders of magnitude in frequency and signal, indicating 
that temporal correlations of 5^(i) decay exponentially 
at B=0. Note this contrasts directly with recent models 
predicting non-exponential dynamics in this regime 0-01 ■ 

Interestingly, the narrow 400 kHz width of the spin 
noise indicates very long hole spin correlation times of 
T^''-400 ns at B=0 {t^ = l/(27rr) for exponential dy- 
namics), greatly exceeds the B=Q hole relaxation 
time measured in InGaAs QDs [l3, [13 ^-^id is consistent 
with the hole relaxation recently inferred at B=Q from 
pump-probe studies of InAs QD ensembles 17| . rj* also 
exceeds the dephasing time of holes localized in narrow 
GaAs quantum wells studied by resonant spin amplifica- 
tion [16|. 

The presence of hole-nuclear coupling becomes plainly 
evident upon applying small transverse fields B^ (see 
Fig. 2). As observed previously 27[, the noise spectrum 
largely shifts to higher frequencies as fluctuations (55^ are 
forced to precess about B^ at the hole Larmor frequency. 
More importantly, however, Fig. 2 also reveals that a 
portion of the zero-frequency spin noise peak remains 
despite application of purely transverse fields. This in- 
dicates that the holes do feel effective nuclear magnetic 
fields in the z direction, because longitudinal fields (real 
or effective) necessarily generate noise at zero frequency. 
In general, spins in an arbitrary magnetic field generate 
two noise peaks: One at high frequency due to trivial 
spin precession, and one centered at zero frequency due 
to longitudinal field components. The former is weak at 
B—Q (for holes) and is strongly broadened due ensemble 
averaging, while the latter is not (for details, see [ll]). It 
is precisely this zero-frequency noise peak that we study, 
as it reveals the intrinsic timescales of (5*^(^)5'' (0)) and 
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FIG. 2. Hole spin noise in purely transverse applied mag- 
netic fields (Bx)- In addition to the expected shift of the hole 
spin noise to the hole Larmor frequency {ul = g\iiBBx/h, 
with g'j^^O.lS), there remains a finite noise component at 
zero frequency. This reveals the presence of the longitudi- 
nal (z) components of the nuclear Overhauser magnetic field, 
Bn. Longitudinal fields, real or effective, necessarily result in 
spin noise at zero frequency [l^]. The integrated noise power 
remains constant. 



the dynamical response function of the hole spin decay. 

To explore the extent to which hyperfine interactions 
limit at 5=0, we apply small longitudinal magnetic 
fields Bz to overwhelm Bj, and suppress hole-nuclear 
coupling. Figures 3(a,b) shows the spectral density of 
hole spin noise as Bz increases to 80 G. The width of 
the noise peak narrows dramatically from 400 kHz to 
less than 40 kHz, indicating that increases over ten- 
fold to nearly 5 /its. More importantly, the detailed line- 
shape of the spin noise evolves away from Lorentzian and 
more closely approaches a 1// power-law decay over the 
measured frequency range, thereby revealing an appar- 
ent crossover from exponential dynamics to much slower 
l/log(t) spin decays in the time domain. Higher fields to 
300 G do not alter the noise lineshape further. 

These data highlight an essential aspect of spin noise 
measurements: the ability to directly reveal detailed 
spectral lineshapes to explore slow and non-trivial decay 
mechajiisms. The data appear to contradict recent the- 
ories predicting slow, inverse-log decays of central 
electron spins coupled to nuclear spin baths at B=Q (we 
see exponential decays at B=Q). Whether these theories 
are fully applicable to holes remains an open question. 
We do observe 1// noise spectra and inverse- log dynam- 
ics emerging in finite (but small) B^ however, suggesting 
at least partial validity of these models. Hole spin decays 
of order 1 jis can arise within models of hole- mediated nu- 
clear CO- flips, but robust Lorentzian noise lineshapes at 
B=Q are not expected (for details, see [11]). One possibil- 
ity is that quadrupolar nuclear interactions and the local 
electric fields in (In,Ga)As QDs could rapidly 'stir' fiuc- 
tuations of Bn at B=Q, accelerating hole relaxation via 
the co-fiip mechanism and leading to exponential decays. 
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FIG. 3. (a) Hole spin noise spectra in purely longitudinal 
applied fields i3z=0, 10, 20, 40, 80 G. The spin noise narrows 
dramatically, indicating over an order-of-magnitude increase 
in hole spin correlation time r]* from 400 ns to ~5 ^s. (b) 
The same spectra on a log-log scale clearly evolve away from a 
Lorentzian lineshape and more closely approach a 1/frequency 
power-law, indicating slow inverse-log spin decays in the time 
domain. The integrated noise power remains constant (as 
expected), (c) as a function of \Bz\. 
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19j or hybridiza- 
have also been proposed, but 
their influence is not explicitly studied here. Note that 
these ensemble studies do not distinguish whether, in this 
regime, every hole in the ensemble exhibits l/log(i) dy- 
namics, or whether there exists a broad distribution of, 
e.g., exponential timescales in the hole ensemble (due 
to varying QD strain and anisotropy), whose sum could 
mimic l/log(t) dynamics. Ongoing efforts are aimed at 
elucidating this difference. 

Although not strictly applicable to power-law dynam- 
ics, if we continue to infer a characteristic timescale via 
the noise half-width over the measured frequency range, 
then Figure 3(c) shows that increases rapidly with in- 
creasing \Bz\, but saturates at 5 /xs when \Bz\ > 100 G. 
These data directly reveal the typical field scale of hole- 
nuclear coupling, ~25 G, or about one-tenth of the 
electron- nuclear coupling in similar QDs [s^, consistent 
with previous studies 15, 2^, 21 1. 

Finally, Figure 4 shows a comprehensive study of 
how temperature and Bz determine and the dynam- 
ical response function. The plots identify three general 
regimes: (i-yellow) At high temperatures T > 15 K, r^' 
falls rapidly, independent of Bz, and noise lineshapes 
are Lorentzian. This suggests straightforward phonon- 
assisted hole spin relaxation mediated by spin-orbit cou- 
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FIG. 4. Identifying different temperature and field regimes of 
hole spin decay. Data are plotted as (a) hole spin correlation 
time versus temperature at various applied fields B^, and 
as (b) Tg versus at various temperatures. Three general 
regimes can be identified (see text). 



pling, and exponential relaxation, (ii-green) Low temper- 
atures below 10 K and in B = 0, single phonon effects are 
suppressed, is limited to 400 ns by hole-nuclear cou- 
pling, and noise spectra remain unexpectedly Lorentzian. 
(iii-blue) Finally, at low temperature using small |i?z| ^ 
80 G to overwhelm Bn, much longer 5 /its hole spin corre- 
lation times are exposed and noise spectra approach 1// 
power-laws, indicative of a crossover to very slow l/log(f) 
decays. 

In summary we have demonstrated that spin noise 
spectroscopy allows unusually detailed studies, at all 
relevant timescales, of dynamic response functions in 
strongly-coupled hole-nuclear spin systems - an inher- 
ently many-body problem that has eluded concise theo- 
retical treatment. Systematic dependencies on tempera- 
ture and magnetic field are revealed, serving as a test-bed 
for theoretical models. The measurement scheme of pas- 
sively detecting the intrinsic spin fluctuations represents 
a kind of 'quantum simulator', which is of great rele- 
vance to other interacting many-body systems of current 
interest including microcavity polariton condensates and 
fractional quantum Hall phenomena. 
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Supplemental Material 



A. Materials and Methods 

Quantum dot structures. InAs/GaAs QDs are grown by molecular-beam epitaxy on (001) GaAs substrates, 
and then thermally annealed at 940 °C [ij. Annealing interdiffuses indium and gallium, resulting in (In,Ga)As/GaAs 
QDs with effective localization volumes of order 2000 nm^, giving ~10^ nuclei within the spatial extent of the resident 
hole's wavefunction. The sample contains 20 layers of QDs, separated by 60 nm GaAs barriers. Each layer contains 
^10^° QDs/cm^. Although not intentionally doped, these QDs are weakly p-type due to background carbon doping; 
we estimate that ^10% of the QDs contain a single resident hole. The inhomogeneously-broadened photoluminescence 
(PL) spectrum of these QD ensembles is typically peaked at ~1.385 eV; see Fig. SI. 

Spin noise spectroscopy. The QD samples are mounted on the cold finger of a small optical cryostat. A linearly- 
polarized continuous-wave probe laser is tuned in wavelength to within the PL spectrum of the QD ensemble and is 
weakly focused through the sample (k||z||n, where n is the sample normal). Stochastic fluctuations of the ensemble 
hole spin projection along the z axis, 6S^{t), impart Faraday rotation fluctuations 69p{t) on the transmitted probe 
laser via the usual optical selection rules for positively-charged trions. Balanced photodiodes detect SOpit), and the 
amplified output voltage SVCt) is continuously digitized and Fourier-transformed to obtain the frequency spectrum of 
the measured noise power [2[. External coils provide longitudinal {B^) and transverse (Bx) applied magnetic fields. 
Background noise densities due to photon shot noise and amplifier noise are eliminated by interleaving and subtracting 
spectra acquired at large Bx (>2000 G), which shifts any spin noise to high frequencies outside the measured range. 
This procedure leaves behind only the noise signals arising from fluctuating hole spins at low fields. Typically the cw 
probe laser power is a few hundred /iW, and it is focused to a rather large (100 /xm) spot on the sample to minimize 
heating and self-pumping of the QDs (see Fig. S2). Crucially, and in comparison with previous work 0, the 
present setup uses low-noise, stabilized probe lasers that now permit accurate and quantitative recovery of the small 
low-frequency spin noise signals that exist in the zero-field limit. 

A consistent measure of the characteristic timescale of the hole spin correlations is obtained from the measured 
half- width F of the spin noise peak that is centered on zero frequency. Specifically, we use = l/(27rr), which is 
precise for Lorentzian noise lineshapes that indicate single-exponential relaxation dynamics. This definition of is 
also used when the noise lineshapes deviate from Lorentzian, even though, strictly speaking, power-law lineshapes 
cannot be characterized by a specific timescale. In this case, F is determined relative to the peak spin noise power 
spectral density that is measured in the lowest frequency bin. 



B. Dependence of integrated spin noise on photon energy of probe laser 

Figure SI shows the inhomogeneously-broadened photoluminescence (PL) spectra of the (In,Ga)As/GaAs QD 
ensemble (solid black lines) under very low excitation conditions by a 1.58 eV (785 nm) laser. This PL arises from 
ground-state recombination of both positively-charged trions (from QDs containing a single resident hole), as well 
as from neutral excitons X'^ (from QDs that are empty). Note that these transition energies are typically very close in 
(In,Ga)As QDs {X'^ being higher in energy by '-^ 2 meV Q) and therefore they overlap in this ensemble PL spectrum 
and cannot be separately resolved. The PL spectrum therefore directly refiects the inhomogeneously-broadened 
distribution of fundamental QD transition energies in the ensemble. 

When performing spin noise spectroscopy of these QDs, the narrow-band, continuous-wave probe laser is tuned 
in energy to within this PL spectrum. Figure SI also shows the frequency-integrated (i.e., total) measured spin 
noise power as a function of the photon energy of the probe laser. The integrated spin noise power provides a 
relative measure of the number of fluctuating spins being measured. Its dependence largely follows the PL spectrum 
with a small blueshift, commensurate with the expected energy difference between X~^ and X'^ transition energies. 
This correspondence indicates that the measured spin noise arises from the resident holes that are trapped in the 
singly-charged subset of the QDs (rather than from spins in, e.g., the buffer or wetting layers or in the bulk of the 
semiconductor wafer). Further, at all probe laser energies where spin noise is detected, the spin noise exhibits the 
same narrow spectral width at zero applied magnetic field (as shown in Figure 1 of the main text), and the measured 
spin noise has the same behavior in transverse and longitudinal fields as shown in the main text (verified for a number 
of different probe energies). 
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FIG . SI. Dependence of the frequency-integrated (i.e. total) spin noise power on the photon energy of the probe laser (points). 
Also plotted is the photoluminescence (PL) spectra of the (In,Ga)As/GaAs QD ensemble (solid line), (a) and (b) show data 
from two different QD ensembles (both annealed at 940 °C). The good correspondence between the PL and total spin noise 
indicates that the spin noise signals arise from resident holes in these QDs, and not from holes or electrons residing in, e.g., 
the wetting or buffer layers of the structure. 

C. Spin noise spectroscopy: a passive measurement of hole spin fluctuations 

The particular implementation of spin noise spectroscopy employed in these experiments to detect fluctuations of 
(the net spin polarization of the resident holes in the QD ensemble) is based on optical Faraday rotation. The 
Faraday rotation angle Op depends on the difference between the indices of refraction for right- and left-circular ly 
polarized light, and vT . In particular, Op{uj) = ^[rt+(w) — n^(aj)], where L is the sample thickness, ui denotes 
energy, and c is the speed of light. 

Before discussing the case of an inhomogeneously-broadened QD ensemble, first consider a spin noise measurement 
of a single homogeneously-broadened optical absorption resonance a(ui) having a Lorentzian line-shape centered at 
ido a-iid half-width 7; namely, a{uj) ^ ; as shown in Figure S2a. Let us say that this absorption resonance 

is spin-dependent as for the case of the optical transition between a resident hole and a positively-charged trion . 
Following the usual optical selection rules, if the resident hole is "spin-up" then it can only absorb a~ circularly- 
polarized light, but if it is in the "spin-down" state then it can only absorb (t~^ circularly-polarized light. 

The dispersive part of this optical transition - that is, the part complementary to the absorption resonances a^{uj) 
- are the indices of refraction n^{uj) ~ (a;-^"^)^°+7^ • ^ function of detuning A = u — cuo from resonance, the 
indices n^{uj) decay much more slowly than the absorption a'^{uj): as A^^ versus A~^, respectively. Therefore, a 
measurement of Faraday rotation can remain sensitive to the spin state of the hole even for large detuning A ^ 7 
where a^{uj) — >■ and the number of photons absorbed by the system becomes vanishingly small. In this regard, the 
measurement of can be considered non-perturbative, and spin noise measurements of this type were performed on 
alkali vapors at large probe laser detunings from the 5 — > P (Dl and D2) lines of potassium and rubidium [s], and 
were also performed on conduction-band electrons in bulk n-GaAs with the probe laser detuned well below the GaAs 
band-edge 0,0. 

Note that 6p and therefore the spin noise measurement is most sensitive to spin fluctuations of this idealized system 
when the probe laser and the absorption resonance are separated by an energy of A ~ 7. The sensitivity falls off 
slowly as A~^ at larger detunings. Even more importantly, note that 6*^- — >■ at small detunings A < 7, and that 
9f = exactly on resonance (A = 0) because there is no difference between n+(a;) and n~(Lu). The spin noise 
measurement is not sensitive to transitions lying at the same energy as the probe laser. 

This brings us to the case of measuring spin noise in an inhomogeneously-broadened QD ensemble. Here we tune 
the laser to directly within the inhomogeneously-broadened distribution of QD transition energies. As shown above in 
Figure SI, this maximizes the measured spin noise signal. Undoubtedly, those QDs in the ensemble that happen to be 
resonant with the probe laser are optically pumped by the probe laser. However as noted in the previous paragraph, 
it is precisely these resonant QDs that do not, to leading order, give any spin noise (moreover, any photogenerated 
electrons, holes and/or trions that are pumped by the linearly-polarized probe laser are not spin-polarized). Rather, 
the probe laser will be primarily sensitive to those resident holes in QDs that have detuned X+ transitions, and 
these QDs are not pumped (see Fig. S2b). In this regard, the Faraday rotation measurement still functions as a 
non-perturbative probe of S^, the net spin polarization of resident holes in the QD ensemble. 
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The fact that 9f{uj = wq) is not sensitive to spins in QD ensembles that are weakly pumped at w = was 
demonstrated in several recent ultrafast studies using independently-tunable pump and probe lasers 0-01 ■ 

Nonetheless, great care is taken to ensure that the probe laser functions only as a passive detector of the resident 
hole spin fluctuations, and does not inadvertently perturb the measurement of S^. Primarily, we wish to ensure that 
the probe laser i) does not heat the QDs which could lead to incorrect spin lifetime measurements, and ii) does not 
inadvertently detect any particles (electrons, holes, excitons or trions) that are photoexcited by the probe laser itself. 
Thus we use low probe laser intensities by using low laser power (~200 /xW) and large spot sizes (~100 /im). Although 
this significantly reduces the measured spin noise signals, it ensures that we operate in a regime where the measured 
spin correlation times {i.e., the inverse width of the spin noise spectra, l/(27rr)) are independent of probe laser 
power. 

In this regime. Figure S2(c) shows that the measured half-width F of the low temperature, zero-field hole spin noise 
spectrum is independent of the probe laser power. Significantly higher probe laser intensity leads to a broadening of 
the spin noise peak, indicating shorter hole spin lifetimes. For the same series of measurements. Figure S2(d) shows 
the total (frequency-integrated) spin noise signal in units of volts of detected signal. As expected, the total spin noise 
signal detects no contribution from photoexcited carriers, since it increases nearly linearly with probe laser power 
(all else being equal, doubling the laser power simply doubles the voltages at the photodetectors). If the probe laser 
were inadvertently measuring spin noise from photogenerated (rather than resident) particles in the resonant QDs, 
the total spin noise would be expected to increase super-linearly (doubling the laser power would not only double the 




FIG . S2. Spin noise spectroscopy using optical Faraday rotation: a non-perturbative measurement of resident hole spin 
fluctuations, (a) Illustrations depicting the a'^ and a~ circularly polarized absorption resonances Q*(a;) and the associated 
indices of refraction (lj) of an idealized homogeneously-broadened (Lorentzian) optical transition {e.g., from a X"*" transition). 
n^{uj) decay much more slowly than a*(tj) as a function of large detuning A. Faraday rotation. Of oc — n~, can therefore 
remain sensitive to the spin state S'^ of the resident hole, even for large A where a* 0. Importantly, note that 9f ~ when 
A = 0. (b) An illustration of the inhomogeneously-broadened distribution of transitions within the QD ensemble. The 
probe laser is most sensitive to hole spin fluctuations in those QDs having X"*" transition energies detuned by about 7 and is 
not sensitive to those QDs that are on resonance (i.e., to those QDs that are unavoidably pumped by the linearly-polarized 
probe laser), (c) The half-width F of the zero-field hole spin noise spectrum measured using different probe laser power: There 
is no perceptible infiuence of the probe laser on the hole spin correlation time, (d) The total (frequency-integrated) spin noise, 
in volts, as a function of probe laser power. The total spin noise scales nearly linearly with the probe laser power, indicating no 
contribution to the spin noise from excited carriers. The dots are measured values, and the solid lines are power-law fittings. 
In (c) and (d) the red and black dots denote two measurements with different spot sizes. 
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voltages at the photodetectors, it would also double - or at least increase - the number of particles being measured, 
leading to a superlinear dependence). 

D. Spin noise measurements of free electrons in n-type bulk Gallium Arsenide 

Figure S3 shows that magnetic fields in the longitudinal (z) direction, either real or eff^ective, necessarily lead to 
some spin noise centered at zero frequency. To demonstrate this, we show spin noise measurements of electron-doped 
bulk GaAs. Free conduction band electrons in n-type bulk GaAs are delocalized and sample a huge number of lattice 
nuclei. Therefore the influence of the fluctuating nuclear spin bath on these free electrons is extremely small. To 
leading order, the only magnetic flelds felt by the electrons are those that are externally applied. In this case, a purely 
transverse applied magnetic field Bx uniformly shifts the spin noise of these electrons out to the Larmor frequency 
i^L — gefJ-sBx/h, and leaves no remnant of spin noise at zero frequency (Figure S3a). This is in marked contrast 
to the spin noise spectra of QD holes shown in Figure 2 of the paper, in which some spin noise clearly remains at 
zero frequency despite application of a purely transverse B^ - this remaining spin noise is due to the longitudinal 
components of the effective nuclear field B„ that is felt by the holes. 




Frequency (MHz) Frequency (MHz) 

FIG . S3. Spin noise from free electrons in bulk n-type GaAs. (a) Electron spin noise in pure transverse fields Bx ~ 0, 2.5, 5.0, 
7.5, and 10 G (black to green). Note the absence of residual noise at Hz, as expected from particles that feel little influence 
of the nuclear spin bath. The probe laser's wavelength is 845.6 nm (well below the GaAs band-edge), its power is 2 mW, and 
the spot size is large (100 /im). This n-GaAs wafer is 350 fj,m thick and is doped at ne — 1.4 x lO^'' cm^"* (it is "sample A" in 
Ref. i). (b) Spin noise in a 60 G canted applied magnetic field (B — |_B|cos^ x+ |i3|sin6' z), showing that longitudinal (z) 
field components lead to spin noise at zero frequency. 

That longitudinal fields generate spin noise at zero frequency is explicitly demonstrated in Figure S3(b), which 
shows electron spin noise spectra in an intentionally tilted applied magnetic field. A 60 G applied field B is rotated 
from the transverse to the longitudinal direction (0 = — ?► 7r/2). As B acquires a longitudinal component, zero 
frequency spin noise grows as sin^0, with a commensurate cos^O suppression of noise signal at the electron Larmor 
frequency. 

E. Spin noise of holes in static nuclear fields 

It is useful to construct a toy model of hole spin noise in the considerably oversimplified limit of static nuclear 
Overhauser fields Bn = {Bn^ , i?n„ :BnJ- The purpose of this exercise is three- fold: 

1) It demonstrates that hole spin precession about Bn generates a broad hole spin noise spectrum at high frequencies 
between 5-100 MHz. The broadness of this high frequency noise is due to the statistical distribution of Bn over the 
QD ensemble, and also to the inhomogeneous distribution of hole g-factors within the QD ensemble. 

2) It demonstrates that the longitudinal Overhauser fields give a delta-function (or at least very narrow) spin 
noise peak at zero frequency, that is not expected to narrow or broaden with applied fields B^ or B^ (in contrast to 
actual experimental observation). This noise peak is not statistically broadened by the ensemble, since each QD gives 
some noise at exactly zero frequency. 

3) It shows that, in zero applied field, the high-frequency precessional noise is strongly suppressed as compared to 
the zero-frequency noise due to the large anisotropy of the hole g-factor. 

In marked contrast with electrons, holes couple very anisotropically to in-plane versus out-of-plane magnetic fields. 
For pure heavy holes, the longitudinal (out-of-plane) g-factor is finite while the transverse (in-plane) g-factor g± is 
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zero. However, hole eigenstates in typical p-type self-assembled III-V quantum dots invariably contain some admixture 
of light hole states in addition to their predominantly heavy-hole character. This leads to a small in-plane g-factor 
g±_ that is of order 0.15 in our QDs, which is about an order of magnitude less than (gy ^ 1). Note also that there 
exists a large inhomogeneous dispersion of these g-factors within the QD ensemble, likely due to differences in QD 
shape and strain. 

Here we assume that the Overhauser field Bn in each quantum dot has components i?„^ , i?„^ , and that are 
each Gaussian-distributed with typical dispersion of ^25 Gauss (taken from experimental data of Figure 3c). In any 
given QD, hole spin precession about the transverse component of generates noise in (the measured quantity) 
at the hole's Larmor precession frequency. 



(where here we have generalized slightly to allow for real applied magnetic fields and/or B^, and we assume no 
in-plane anisotropy of g± for simplicity) . 

In each QD, ujl occurs at a different frequency depending on the magnitude and direction of Bn in that QD, and 
also depending on g± and g\\ in that QD. Using the typical values stated above, lol can range from a few MHz out 
to ^100 MHz. Averaging over many QDs, this precessional noise generates a very broad spectrum, weakly peaked at 
about 10 MHz and spanning the range from 5-100 MHz. It is this precessional noise that represents the trivial ensemble 
dephasing of an ensemble of hole spins that are all initially oriented at t = 0, as in a pump-probe measurement. 
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FIG . S4. A toy model of hole spin noise in a QD, assuming a static nuclear field Bn. In each QD, hole spin precession about 
Bn induces a high-frequency noise peak at the hole Larmor frequency ujl- In addition, longitudinal (z) components of Bn 
generate a noise peak at zero frequency (for the analogous case of electrons in arbitrary applied fields, see Fig. S3 above). For 
the purposes of this toy model, these peaks are represented by delta functions. The total noise power contained within the 
high-frequency peak is generally much smaller than that within the zero-frequency peak, due to the large anisotropy of the hole 
g-factor. Statistical averaging over the QD ensemble smears the high-frequency precessional noise over a large frequency range. 
On the other hand, the zero- frequency noise peak is contributed to by every QD. In our noise experiments, it is precisely the 
width and lineshape of this zero-frequency peak that we seek to measure, as these parameters reveal the intrinsic spin correlation 
timescales and decay mechanisms of the hole spins. 



However, in each and every QD, the longitudinal component of Bn {i.e., Bn,), generates a finite time-averaged 
value of the hole's spin projection, (5^(t)), that does not decay in time. In a noise measurement, and within the 
limitations of this toy model, this simply gives a delta function at zero frequency, and applied magnetic fields do not 
alter the lineshape of this peak. Of course in reality, {S^{t)) does decay because Bn is not static, and the associated 
noise peak is not a delta function - it is precisely the width and lineshape of this zero-frequency peak that reveals the 
long-time decay mechanisms of the central hole spins that we are interested in. 

Within this simple model, it can be shown that the total (frequency-integrated) power of the noise peak at ujl is 



gi(B, + B„j2+.gi(B„j2 



while the integrated power of the noise peak at zero frequency is 



(2) 



° gl{B.:, + BnJ^+gl{B,J^+gl{B^ + B^J^- 
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Therefore in the absence of any appUed fields {B^ = = 0), the noise power at wl is typicaUy suppressed as 
compared to the zero- frequency noise by ~ {g±/g\\)'^ (except when _B„^ is very smaU). This, together with the above- 
noted fact that ensemble averaging smears out the high-frequency precessional noise over a very broad bandwidth, 
may explain why we do not observe any clear sign of precessional noise occurring at high frequencies in the absence 
of applied fields. Only the much larger zero-frequency noise is apparent. However, as we apply large magnetic fields 
in the transverse direction {Bx), all the hole noise power is expected to shift to the precessional noise component at 
wl, in agreement with our experimental data (Figure 2 of the main paper). Further, as we apply large fields in the 
longitudinal direction (B^), the total noise power contained within the zero- frequency component will increase only 
very slightly (all the noise power is essentially already contained in this peak), again in agreement with experimental 
observation. 

Finally, note that the above expressions also hold for electrons, in which case g± ^ and the high-frequency and 
zero-frequency noise peaks are expected to have comparable integrated power. 



F. Contribution to the hole spin noise from hyperfine-mediated "co-flips" of nuclear spins 



Previous theoretical studies of nuclear spin "co-flips" (a flip-flop of two distant nuclei mediated by virtual hole- 
nuclear spin flips) revealed non-exponential decay of the central spin due to non-Markovian evolution and strongly- 
correlated dynamics 13, 11|. 



The goal of this section is to show that the time scale of ^1 /iS, which roughly corresponds to the inverse width 
of the spin noise power in zero magnetic field, is in agreement with a simple theory of how nuclear co-flips lead to 
central spin relaxation. We will also highlight differences between this theory and the observed behavior of noise 
power spectrum at longer time scales. 

Consider the effective Hamiltonian of a single hole spin interacting with N ^ 10^ nuclear spins, in the presence of 
an applied out-of-plane (longitudinal) magnetic field B^- We have 

H = LuJ, + V, (4) 



V = '^Yj_Si.xSx, = .g||S^ + ^7^|S,^, (5) 

i i 

where S stands for the hole spin operator, s; is the i-th nuclear spin operator, 7|| and j± are the coupling strengths, 
out-of-plane and in-plane respectively, between the central spin and nuclear spins. We assume, for simplicity, that 
all spins are 1/2, and 7|| and 7^ are in the same ratio as the ratio of longitudinal and transverse hole g-factors, 
and g±. Thus, 7|| is about ten times greater than 7^. Note that in general 7^ is complex, reflecting the fact that 
transverse coupling involves both x and y spin components, but as this does not influence the following discussion we 
treat these couplings as real positive parameters. 

The large value of N and the coupling anisotropy enable a perturbative approach [10], in which the zcroth order 
wave function of a typical state can be approximated to be the eigenstate of the total spin operator along z-axis, for 
example. 



1*) « Ivf-") = |^,;...;t...t), (6) 

where "ff or JJ- indicates the state of the hole central spin (up or down along z-axis) and other arrows show states of 
the nuclear spins. We will call the expectation value of operator uiz in the state j^*) the bias, and we will denote it by 
U!z- We will show that this part of the hyperfine coupling resists the flip-flop transitions. 

Diffusive dynamics of the Overhauser fleld Bn are possible because the couplings of central spin to different nuclear 
spins are not equal so that exchanging direction of a pair of nuclear spins changes the bias value by the amount 
~ 7||. According to the Hamiltonian (j4|), such pair- wise co- flips of nuclear spins happen due to transitions through 
the virtual states with a flipped central spin. The accumulation of pair-wise nuclear spin-flips leads to diffusion in the 
space of Overhauser fleld values. The expectation value of the central spin follows the direction of the Overhauser 

fleld. When the bias changes across the region with B^ + Bn^ ^ {ia-Ii\\)\J (Bn^), the variation of this direction 
becomes substantial, corresponding to the relaxation of the central spin expectation value. 
There is no direct coupling between the state ([6]) and the state 



(7) 
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where nieans that nuclear spins i and j change their directions in comparison to the state j^f). Therefore, to 

estimate the rate of nuclear spin co-flips we use first order perturbation theory to incorporate transitions to the virtual 
states in our wave functions: 

\^)^c(^\^")-J2^jV^)^ (8) 

where 

= (9) 

and where {i} means that i-th nuclear spin is flipped from up to down state, and C ensures a normalization. This 
definition of a typical state j^") includes fast transition processes to the virtual states \V^). One can see that V already 
directly couples states IvP) and 

Further perturbation expansion breaks down because there are many states l^^ij) with almost the same energy as 
the initial state \^). The presence of such quasi-degenerate states leads to incoherent transitions from jvl/) into one 
of the states The rate 1/r of such transitions can be estimated by Fermi's Golden Rule as follows: There are 

N'^ /A pairs of distinct states (again, we are considering spin 1/2 nuclei for convenience), distributed around the 

mean value with interval of energies i57||, where i57|| is the width of the distribution of around their mean value. 
Thus the rate, t^^, of an incoherent transition from the initial state to one of the states that differ from 

1 5') by CO- flips of two nuclear spins is given by 



where iV^ / (4i57| | ) is the density of states represented by all possible pairs and 7^ /4:UJz is the strength of a typical 

couphng between states l^^ij) and the state j^*), i.e. {^ij\V\'^). 

We use 7|| ~ Mflo/Q, where according to various estimates, M ^ 3-13 fieV 12, l^] is the hyperfine couphng per 



nuclear spin per unit probability of the hole being in the unit cell that contains this nuclear spin. The factor flQ/fl 
is the ratio of volumes of the unit cell Oq ^ (0.57 nm)'^ to the volume of the quantum dot, ^ 2000 nm"^. Hence 
no/n - 10-* and 711 ~ 1 • 10^ s'^. 

The fastest spin relaxation happens in states that have initial bias < ^^\/N /2 because such states experience 
comparable couplings along longitudinal (along z-axis) and transverse axes. In such states, changes in total field along 
z-axis (also the bias lOz) leads to a substantial change of the orientation of the total field that acts on the hole's spin. 
If we consider a state with initial lOz close to lo^ ^ ^^\/N /2 we find from (jlOp that the typical time of an incoherent 
transition between states that differ by a co-flip energy is r ~ 10~^ s. 

This time r should not be confused with a hole's spin relaxation time because a single co-flip does not lead to 
a substantial change of the total hyperfine field. Moreover, such random transitions can lead both to an increase or 
a decrease of uJz- However, accumulation of co-fiip transitions leads to a diffusion in the space of bias values. Hence 
the latter changes with time according to the diffusion law: 



where the coefficient 7| | reflects the fact that a single co-fiip corresponds to a change of the bias of order 7[ [ . The 
hole's spin rotation angle becomes of order unity when the change of the bias becomes comparable to the initial bias, 
7x\/]V /2. This gives us an estimate of the hole's spin relaxation time: 

The experimentally obtained value, ^ 400 ns, is in very reasonable agreement with this estimate. Here we note, 
however, that the region of bias values with Uz < 7±ViV' /2, accounts for only about 10% of statistically possible states 
of the nuclear spin bath. When ujz exceeds the size of this region, incoherent flip-fiop transition rate is suppressed, as 
it ca n be seen from the fact that uJz enters as a second power in the denominator in Equation ([T0| . Thus, for values 
of (uJz) ~ 7||'\//V /2, the fiip-flop transition rate is suppressed by a factor 100, in disagreement with the assumption 
of a single relaxation time for all quantum dots, which is expected from the Lorentzian shape of the power spectrum. 
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Dipole interactions between nuclear and central spins allow additional coupling terms in the interaction Hamiltonian, 
such as ~ SzSix that could contribute to the nuclear spin dynamics, but such terms were estimated to be negligibly 
small for heavy- hole states in quantum dots 12| . It is also possible that strains and nonuniform doping introduce high 
gradients of local electric fields that couple to quadrupole moment of nuclear spins. Quadrupole nuclear interactions 



with nonuniform electric field often broaden the solid state NMR-lineshapes of nuclei to the MHz range [IJ]. Such 
interactions do not directly involve the central spin but, if sufficiently strong, they can induce fast intrinsic dynamics 
of nuclear spins and stir fluctuations of the Overhauser field across its typical values. It is possible that, supplemented 
by fast incoherent co-flip processes near the zero bias values, such fluctuations lead to the relaxation of the central 
spin at a fraction of a microsecond. 

Finally, we discuss consequences for the case of nonzero applied fields B^- According to our theory, the window 
of Overhauser field values that allow fast spin relaxation due to co- flips corresponds to \Bz + -BnJ < "y±\/N- The 
average number of quantum dots that happen to have the hyperfine field inside this narrow window at same moment 
of time follows the Gaussian distribution of the hyperfine field |B„|. Consequently, this number is suppressed when 

l-Bzl > '\J (Bn^), while the relative contribution to power spectrum from the states having the bias far from the fast 
relaxation window increases with Bz- Hence, the width of the noise power spectrum should quickly decrease with Bz 



when \Bz\ > \ (B„ ), in agreement with the experimental data. 
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